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Abstract 

We employ QCD sum rules to calculate the J/ip form factors and coupling constants with 
charmed mesons, by studying three-point correlation functions. In particular we consider the 
J/ipD*D and J/ipDD vertices. We determine the momentum dependence of the form factors for 
kinematical conditions where J/ip, D or D* are off-shell. Extrapolating to the pole of each one of 
the so obtained form factors, we determine the coupling constants. For both couplings our results 
(within the errors) are compatible with estimates based on constituent quark model. 
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I. INTRODUCTION 



Since their begining in 1998, effective Lagrangian theories for charmed mesons were con- 
structed to give us a better understanding of the interaction between J / ip and light mesons 
in nuclear matter. Although the original motivation is still very strong, specially in view of 
the vigorous ongoing RHIC program, this new sub-area of hadron physics introduced new 
questions which are now in debate. Some of them are: What are the coupling constants 
involving heavy mesons? Which are the form factors appearing in these vertices? To what 
extent are we allowed to use SU(4) symmetry? Moreover, interactions of charmed mesons 
offer a testing ground for the predictions of heavy quark effective theory (HQET) 1]. As 
pointed out in ref. Q], a consequence of the spin symmetry of HQET is that 

_ g^DD n \ 

94>dd* — • UJ 

m D 

We will in this work calculate both couplings independently and will have the opportunity 
to check the relation above. 

Charmed hadrons are composites of the underlying quarks whose effective fields describe 
point-like physics only when all the interacting particles are on mass-shell. When at least one 
of the particles in a vertex is off-shell, the finite size effects of hadrons become important. 
Therefore, the knowledge of form factors in hadronic vertices is of crucial importance to 
estimate any amplitude using hadronic degrees of freedom. This work is devoted to the 
study of the J/ipD*D form factor, which is important, for instance, in the evaluation of the 
dissociation cross section of J / tp by pions and p mesons using effective Lagrangians DDD 

The J/ipD*D coupling has been studied by some authors using different approaches: 
vector meson dominance model plus relativistic potential model |3| and constituent quark 
models Q . Unfortunately, the numerical results from these calculations may differ by almost 
a factor two. The relevance of this difference can not be underestimated since the cross 
sections are proportional to the square of the coupling constants. In ref. it was shown 
that the use of different coupling constants and form factors can lead to cross sections that 
differ by more than one order of magnitude, and that can even have a different behavior as 
a function of yfs. 

In previous works we have used the QCD sum rules (QCDSR) to study the D*Dtt 
DDp [8] and J/ipDD [9] form factors, considering two different mesons off mass-shell. In 
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these works the QCDSR results for the form factors were parametrized by analytical forms 
such that the respective extrapolations to the off-shell meson poles provided consistent 
values for the corresponding coupling constant. In this work we use the QCDSR approach 
to evaluate the J/ipD*D form factors and extend the procedure described above to estimate 
the J/ipD*D coupling constant considering each one of these mesons off-shell, calculating 
the associated form factor, and then extrapolating the three obtained form factors to the 
corresponding on-shell masses. We will apply this method to both J/tpDD and J/ipD*D 
vertices. 



II. THE QCD SUM RULE APPROACH 

The three-point function associated with the hadronic vertex J/tpDD*, where J/ip and 
D are the incoming and outgoing external mesons with momentum p and p' respectively 
and D* is the off-shell meson, with momentum q = p' — p is given by 

r»AP,P) = J d 4 xd 4 y e**e^*(0\T{j D {x)j?\ y )jf(0)}\0) , (2) 

where the interpolating fields are = cry^c, j®* = qj v c and jo = iqi^c with q and c being 
a light quark and the charm quark fields. 



The fundamental assumption of the QCD sum rule approach [10|,C| is the principle of 
duality, i.e., it is assumed that there is an interval over which a correlation function may be 
equivalently described at the quark level and at the hadron level. Therefore, the procedure 
of the QCD sum rule technique is the following: on one hand we calculate the correlator 
in Eq.([2|) at the quark level in terms of quark and gluon fields. On the other hand, the 
correlator is calculated at the hadronic level introducing hadron characteristics, giving sum 
rules from which a hadronic quantity can be estimated. 

The hadronic side of the correlation function, r^^p'), is obtained by the consideration 
of J/if>, D* and D state contributions to the matrix element in Eq. (J2J: 

r<^ W) = Hgfiw w , {q2 gff^gl ml) + h. , , (3, 



where h. r. means higher resonances. 

The meson decay constants appearing in Eq. (j^J) are defined by the vacuum to meson 
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transition amplitudes: 



(0\j^\D) = ^R, (4) 



and 

(V(p,e)\jl\0) = m v f V f? at (5) 

for the vector mesons V — J ftp and V = D*. The form factor we want to estimate is defined 
through the vertex function for an off-shell D* meson: 

<V>(p, \)\D(p')D*(q, A')) = 4^1(?V^4(P)<(<?H% • (6) 

The contribution of higher resonances and continuum in Eq. © will be taken into account 
as usual in the standard form of ref. 

The QCD side, or theoretical side, of the correlation function is evaluated by performing 
Wilson's operator product expansion (OPE) of the operator in Eq. (j2J). Writing V ^ u in terms 
of the invariant amplitude: 

rV(p,p') = A(p 2 ,p' 2 , ? 2 )e Q ^ a / , (7) 

we can write a double dispersion relation for A, over the virtualities p 2 and p' 2 holding q 2 
fixed: 

A(pVV) = / dsdu P ^ q2) , (8) 

Att J (s — p z ){u — p ) 

where p(s,u } q 2 ) equals the double discontinuity of the amplitude A(p 2 , p' 2 , q 2 ) on the cuts 

Am 2 < s < oo, m 2 < u < oo, We consider diagrams up to dimension three which include 

the perturbative diagram and the quark condensate. 

To improve the matching between the two sides of the sum rules, we perform a double 

Borel transformation in both variables P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — > M' 2 . We get 

the following sum rule: 

C ^BEL - e -TF* e ~*& = / / c?sdw p {D *\u, s,t)e~^e'^Q(u max - u) 

(t-m 2 D ») AlT 2 Jim? Ju mm 



m c (g 2 G 2 ) 



+ "TV L n G (t,M 2 ,M 12 ), (9) 



2vr 

where C = !! 7 ^ £ m D */ Z )*m v ,/ v „ t = g 2 , 



3m r / _ s\2 



P^Ku,s,t) = ^[l + -y), (10) 



with A = (u + s — t) 2 — 4us, X2 = u + t — s + 2m 2 and 

1 



max 



-st + m 2 c (s + 2t) ± Js(s - im 2 )(t - m 



(11) 



The last term in Eq. Q gives the gluon condensate contribution. The full expression for 
the gluon condensate contribution is given in Appendix A. 

In Eq. (jHI) we have transferred to the QCD side the higher resonances contributions 
through the introduction of the continuum thresholds So and Uq. The sum rules for the form 
factors considering the mesons J/ip and D as off-shell mesons can be obtained in a similar 
way. In the case that the D meson is off-shell we get basically the same sum rule, the only 
difference being the exchange between m 2 D and m D * in the left hand side of Eq. Q. The 
spectral density is exactly the same Q]. In the case that J/ip is off-shell we get 

(j/i/>) //\ 22 -, 
c g^W o^ = 1 1 1 dsdup (J^ {UjS}t)e -^ e ~^Q( Umax -u), (12) 

(t — m±) 4^ Jm 2 Jur, 



ni- 1 n 



with 



p {JM {u, s,t) = ^l \{u - s) 2 -t(u + s- 2m 2 c )} - Art 2 < qq > 5{s - m 2 c )5{u - m 2 c ), (13) 



A3/2 

and 



max _ 

ruin ~ 2m 2 



st + m 2 c (2s + t)± Jt(t - 4m 2 c )(s - m 2 ) 2 



(14) 



III. FORM FACTORS IN THE J/ip VERTICES 
A. J/il)DD* 

The parameter values used in all calculations are: m D = 1.87 GeV, m D * = 2.01 GeV, 
= 3.1 GeV, fj/$ = (405±15) MeV :14]. For fp and fo* we use the values evaluated in the 
two-point QCD sum rules under the same kind of approximations [15]: fo = (170±10) MeV, 
fo* = (240 ± 20) MeV. For the charm quark mass, quark condensate and gluon condensate 
we use the values normally used in QCD sum rules calculations llj"i c = 1.3 GeV, 
(qq) = -(0.23) 3 GeV 3 , (g 2 G 2 ) = 0.5 GeV 4 For the continuum thresholds we take s = 
{m^ + A s ) 2 and uq = (mom*) + A u ) 2 for the sum rules when D(D*) is off-shell, and 
So = (mo* + A s ) 2 and uq = (mo + A u ) 2 for the sum rule when J/ip is off-shell. We take 
A S = A U = (0.5 ±0.1) GeV. 



FIG. 1: M 2 and M' 2 dependence of gfj? D *{Q 2 = 0). 



We first discuss the J/i/jDD* form factor with an off-shell D* meson. Fixing Q 2 = —q 2 = 
and A s = A u = 0.5 GeV we show in Fig. 1 the Borel dependence of the form factor 
9tfDD*(®)- We see that we get a very good stability for the form factor as a function of 
the two independent Borel parameters in the considered Borel regions. The same kind of 
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FIG. 2: M' 2 dependence of perturbative contribution (dashed line) and the gluon condensate 
contribution (dotted line) to 9^^ DD * a t Q 2 = 2.0 GeV 2 . The solid line gives the final resit for the 
form factor. 
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stability is obtained for other values of Q 2 and for the other two form factors 



In Fig. 2 we show the perturbative (dashed line) and the gluon condensate (dotted line) 
contributions to the form factor g^fjJ DD * (Q 2 ) at Q 2 = 2 GeV 2 as a function of the Borel mass 
M at a fixed ratio M' 2 /M 2 = m 2 D lm\. We see that the gluon condensate contribution is 
negligible, as compared with the perturbative contribution. The same kind of behaviour is 
obtained for other values of Q 2 . This is a very interesting result since it might be indicating 
a convergence of the OPE. 

In Fig. 3 we show the behavior of the form factors gj/^DD*(Q 2 ) at Q 2 = 2 GeV 2 as a 

2 f-D*) 2 9 

function of the Borel mass M' . The solid line gives gj/^oD* a ^ a fi xe d ratio M' /M = 
m z)/ m ^- The dashed line gives g[JJ^ DD * at a fixed ratio M' 2 /M 2 = m^/m^,, and the dotted 
line gives g[//^ D * at a fixed ratio M' 2 /M 2 = m 2 D /m 2 D ». 

4.5 i ■ 1 ■ 1 ■ 
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FIG. 3: M' 2 dependence of 9 y ^j^ DD * (solid line), 9^^ DD * (dashed line) and 9 ( j/^d D * (dotted line) 
for Q 2 = 2.0 GeV 2 . 



AD*) 



We can see that the QCDSR results for g j/^ DD * and g ( u^ DD * are very stable in the interval 
3 < M' 2 < 6 GeV 2 . In the case of g^j/^ijD* ^ ne stability is not as good as for the other form 
factors, but it is still acceptable. 

Fixing M 2 and M' 2 at the values of the incoming and outgoing meson masses we show, 
in Fig. 4, the momentum dependence of the QCDSR results for the form factors g^p*, 
Q^bDD* an d gipDD* through the circles, squares and triangles respectively. Since the present 
approach cannot be used at Q 2 « 0, in order to extract the g^DD* coupling from the form 



(D) 




FIG. 4: Momentum dependence of the J/ipDD* form factors. The dotted, dashed and solid lines 
give the parametrization of the QCDSR results (triangles, squares and circles) through Eqs. (|15[). 
(fTB|) and (fTT)) respectively. 

factors we need to extrapolate the curves to the mass of the off-shell meson, shown as open 
circles in Fig. 4. 

In order to do this extrapolation we fit the QCDSR results with an analytical expression. 
We tried to fit our results to a monopole form, since this is very often used for form factors, 
but the fit was only good for glj^*- For g^ D * and g^jj* we obtained good fits using a 
Gaussian form. We get: 

JJM ( q2) _ 200 _ 4 2 57 ~ m ' (15) 

g$UO?) = W (16) 
gfUQ")^e-^. (17) 

These fits are also shown in Fig. 4 through the dotted, dashed and solid lines respectively. 
From Fig. 4 we see that all three form factors lead to compatible values for the coupling 
constant when the form factors are extrapolated to the off-shell meson masses (open circles 
in Fig. 4). 

All results showed above were obtained using A s = A u = 0.5 GeV. In Fig. 5 we use 
the form factor g^DD*(Q 2 ) to ilustrate the dependence of our results with the continuum 
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thresholds. 
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FIG. 5: Continuum threshold dependence of the form factor g^£,*{Q 2 )- The dashed solid and 
dotted lines give the parametrization of the QCDSR results for A s = 0.4 GeV, 0.5 GeV and 
0.6 GeV respectively. The lower, intermediate and upper set of curves show the results for Z\ w — 
0.4 GeV, 0.5 GeV and 0.6 GeV respectively. 

In Fig. 5 there are three sets of three curves. They show the extrapolation of the QCD sum 
rule results, through parametrizations similar to the one given in Eq. (|17j) . to all kinematical 
region. The lower, intermediate and upper sets were obtained using A u = 0.4 GeV, 0.5 GeV 
and 0.6 GeV respectively. The dashed solid and dotted lines in each set were obtained using 
A s = 0.4 GeV, 0.5 GeV and 0.6 GeV respectively. We see that the dispersion in the region 
< Q 2 < 4.5 GeV 2 , where we have the QCDSR points, does not lead to a bigger dispersion 
at Q 2 = —mjj, where the coupling constant is extracted. Therefore, the extrapolation 
procedure used here does not act as a lever effect for the uncertainties. 

Considering the uncertainties in the continuum threshold and the difference in the values 
of the coupling constants extracted when the D, D* or J/ip mesons are off-shell, our result 
for the J/^DD* coupling constant is: 

9ipDD* = (4.0 ± 0.6)GeV- 1 . (18) 

From the parametrizations in Eqs. (|15|). ()16|) and (fTTj) we can also extract the cutoff 
parameter, A, associated with the form factors. The general expression for the Gaussian 
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parametrization is: Aexp [— (Q 2 + B) 2 /A 4 ], which gives A ~ 4.5 GeV when the off-shell 
meson is D or D*. For the monopole parametrization the general expression is: g [(A 2 — 
m 2 )/(A 2 + Q 2 )]. Therefore, for an off-shell J/ip we get A ~ 7.5 GeV. It is very interesting 
to notice that the values of the cutoffs obtained here follow the same trend as observed in 
refs. : the value of the cutoff is directly associated with the mass of the off-shell meson 

probing the vertex. The form factor is harder if the off-shell meson is heavier, implying that 
the size of the vertex depends on the exchanged meson. Therefore, a heavy meson (like J/i/j) 
will see the vertex approximately as point like, whereas a light meson will see its extension. 

In Table I we show the results obtained for the same coupling constant using different 
approaches. While our result is compatible with the coupling obtained using the constituent 
quark meson model [2], it is half of the value obtained with the vector meson dominance 
(VMD) model plus relativistic potential model jjj. It is important to mention that in the 
case of the VDM model, the coupling is extracted for an off-shell J/ip with Q 2 = 0. However, 
since our g^DD*{Q 2 ) form factor depends weakly on Q 2 , the value extracted at Q 2 = is 
still very close to the value in Eq. (fTKj). 



B. J/^DD 

In ref. |9[ we have studied the J/ipDD vertex, and we have calculated the form factors 
considering two cases: i) one D and ii) the J/ip as off-shell mesons. We obtained 

n ( J M(n 2 ) - 306 - ti 7 63 Z m ^ h<rt 
^ W, "63 + Q 2 " 5, 63 + Q 2 ' 1 Jj 

Cn(Q 2 ) = lSe-^, (20) 

which lead to the coupling 

gi>DD = 5.8±0.9. (21) 

The numerical values presented here are a little different from the values in ref. 0, due to 
the factj;hat there we made a mistake in the value of fj/^, and have used fj/^ = 270 MeV 



in ref. 
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IV. CONCLUSIONS 



We have used the method of QCD sum rules to compute form factors and coupling 
constants in the J/ipDD* and J/ipDD vertices. We have first analyzed the J/ipDD* vertex 
and have considered three cases: i) off-shell D, ii) off-shell D* and iii) off-shell J/ip. In 
the three cases we have fitted the QCDSR results with analitycal forms and extracted the 
coupling constant. Our results for the coupling show once more that this method is robust, 
yielding numbers which are approximately the same regardless of which particle we choose 
to be off-shell and depending weakly on the choice of the continuum threshold. As for the 
form factors, we obtain a harder form factor when the off-shell particle is heavier. The same 
comments can be made for the J/ipDD vertex. 

In this work we have not considered a s corrections, which might be not negligible in the 
case of heavy quarks. However, as shown for instance in ref. in the calculation of B*Btt 
and D*Dtt coupling constants, a s corrections are more important in the beauty case than 
in the charm case. In their caculation the inclusion of a s corrections have changed gD*D-K 
by about 10%, which is smaler than the theoretical uncertainties. Therefore, we postpone 
the evaluation of the a s corrections to a future work. 

As a closing remark we come back to Eq. (|21|) from where we get g^Doj^D = (3.1 ± 
0.5) GeV" 1 , which is in agreement (considering the uncertainties) with our result for g^DD* 
in Eq. (fTHj). Therefore, our QCDSR results for the couplings g^DD and g^DD* obey the 
HQET relation given in Eq. 

In Table I we present our final results and compare them with other calculations. 



coupling 


this work 


ref. [2] 


ref. [2] 


g^DD* (GeV)" 1 

g^pDD 


4.0 ± 0.6 
5.8 ± 0.8 


8.0 ± 0.6 
7.7 


4.05 ± 0.25 
8.0 ± 0.5 



TABLE I: Values of the coupling constants evaluated using different approaches. 

We see that while our result for g±r>D* is compatible with the coupling obtained using 
the constituent quark meson model |2j, this is not the case for g^oo- The values of the 
couplings obtained with the vector meson dominance model plus relativistic potential model 
(jil are bigger than ours for both couplings. The origin of these discrepancies deserves further 
investigation. 
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APPENDIX A: GLUON CONDENSATE CONTRIBUTION 
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FIG. 6: Gluon condensate diagrams contributing to the correlation function in Eq. 



In Fig. 8 we show the gluon condensate diagrams contributing to the correlation func- 
tion in Eq. (j2J). We give bellow the expressions for all diagrams contribution in terms of 
II G (t, M 2 ,M' 2 ) appearing in Eq. ©. 

1 



n { e\t,M 2 ,M' 2 ) 



A//"' 6 f OO , n 2 a(M 2 + a M' ) t (aM ,2 -l) 

" 2 7 »,r/2 



M 2 (M 2 + M' 2 ) 2 A/M' 2 



G?tt e M 2 (aM' 2 -\) g M 2 +M' 



+ 



Tm 2 («m' 2 - 1) 



(Al) 



GV ' ' 7 6M 2 (M 2 + M' z ) Ji/M' 2 

J ( — M' 2 + 2aM' 4 - M 2 (3 - 5aM' 2 + a 2 M' 4 ) - m 2 a(M 2 + M /z (2aM' z + 

aM' 2 - 1 V c 



r/2\ 



M 



/2 



+ 



m^(M 2 + aM l4 ) t(oM /2 -l) 
j^q, g M ! (aM' J -l) g M 2 + M' 2 

2 



atM' 2 + 



/2/ 



r/2 



, w2 m 2 a 2 M' z (M 2 (M 2 + 2M' Z ) - aM'°(aM' A - 2)) , r2 , . 4 
1 + aM' —2— - — ^ i — — ^ l -L + M + 3aM' + 



M 2 



aM' 



+ 



tM ,2 {aM 2 + 1)(M 2 + aM' 4 ) M' 2 (M 2 + aM' 4 )(2aM' 2 + aM 2 - 1)) 



(M 2 + M' 2 ) 2 



(M 2 + M' 2 )(aM' 2 



(A2) 
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n^(t, M 2 , M' z ) = 0, since it is proportional to m q , 



r/2 



(A3) 



ng 1} (t, m 2 ,m' 2 ) 



/2a 



M' 2 



6M 2 Ji/m' 2 



m 2 a (M 2 + aM' 4 ) t( Q M ,2 -l) 
(^Q/ g M 2 (aM' 2 -l) g M 2 + M' 2 



ft 



(aM" - 1) L 



m\aM' 2 2{M 2 + aM'*) \ (m 2 c a 2 M'\M 2 - aM'* + 2M' 2 ) tM''\aM 2 + 1) 



,2^,2 /\/f/2 



/■I 



r/2 a 



r/2/ 



+ 



aM' 2 - 1 
/8(M 2 + M' 2 ) 
2(oM' 2 - 1) V aM' 2 -l 



M 2 (aM' 2 - l) 2 



+ 



(M 2 + M' 2 ) 2 



/2 



am 2 M' 2 



(A4) 



njf } (£,m 2 ,m 



/2\ 



M 



/i 



12M 2 7i/m' 2 



m^atA^ + aM' 4 ) t(aM /2 -l) 
^q, g Af 2 ( Q M' 2 -l) g M 2 + M' 2 . 



a z 



a 



M' 2 - l) 2 



'm 2 c a 2 M' l [M 2 + M /2 )(M 2 + 2M / 2 - aM' 4 ) tM' 2 (aM 2 + 1) | 



M 2 («M' 2 - l) 2 



(M 2 + M' 2 ) 



(M 2 + 2M' 2 + aM' 2 M 2 ) \ ( ' m 2 c a 2 M' 2 (M 2 + 2M' 2 - aM' 4 ) tM' 2 (aM 2 + 1) 



+ 



a 

m 2 c a 2 M' 2 (M 2 + M' 2 )(aM' e - aM' 4 -3aM' 2 M 2 -AM' 2 - M 2 ) { m 2 c a 2 M 



M' 2 - 1 



M 2 (aM' 2 - 1) 



+ 



r/6 



/2; 



r/2 



(M 2 + M' 2 ) 2 

,2^,2^/4 



(M 2 - aM' 4 + 2M 



M 2 (aM' 2 - l) 3 
tM 2 M' 2 (aM 2 + 1) 



+ 



M 2 



+ 



+ 



(aM' 2 -l) 2 / (M 2 + M /2 ) 2 
M /2 ((aM 2 + 2)(2aM' 2 + 1) + 3aM 2 )' 
(aM' 2 - l) 2 ^ 



2taM'\aM 2 + 1) 
(M 2 + M' 2 )(aM' 2 - 1) 



+ 



+ 



(A5) 



n 



f\t,M 2 ,M 



/2a 



m 2 a 2 M' 4 

1 + 



M' 



/2 



3 Jl/M' 2 
r2/„.)i^/2 



(^Q/ g M 2 (aM' 2 -l) g M 2 + M' 2 



a 



(aM' 2 - 1) 



tM 2 (aM' z - 1) 
M 2 (aM' 2 - 1) ~ (M 2 + M' 2 ) 2 



(A6) 



ng 4) (t,M 2 ,M ,z ) 



/2a 



1 



hi 



6M 2 Ji/m' 2 



m 2 a(M 2 + aM /4 ) t(aM /2 -l) 
(lOl g M 2 (qM' 2 -1) g M 2 + M' 2 . 



M' 2 - 1) . 



a 



+ 



my(M 2 + aM' 4 ) taM' 2 (3M 2 + M' 2 — aM 2 M' 2 + aM' 4 ) m 2 c ta 2 M' 2 

+ (M 2 + M' 2 ) 2 + 

M 4 (-6 + 14aM' 2 - 9a 2 M' 4 ) + M 2 M /2 (-2 + 3aM' 2 + 



2(aM' 2 - 1) 
1 



(M 2 + M' 2 )(aM' 2 - 1) 
3a 2 M' 4 ) + M ,4 (-l - a 

+ 2M 4 (-1 + 3aM' 2 ) + M 2 M ,2 (1 + 6aM' 2 + a 2 M' 4 ))) 



- 3a 2 M' 4 ) + M ,4 (-l - aM' 2 + a 2 M' 4 ) + +m 2 a(M' 4 (l + 4aM' 2 - a 2 M' 4 ) + 



(AT) 
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ng } (t, m 2 , m' 2 ) = ng x) + ng 2) + ng 3) + ng 4) ( as) 



U^\t,M 2 ,M ,z ) = 



/2\ 



i2 



QM 2 (M 2 + M 1 ) Ji/M' 



mgo(M 2 +«M' 4 ) t( a M /2 -l) /\/f2 _|_ q,^' 4 



(/q/ g M 2 (nM' 2 -l) g M 2 +M' 2 



«M' 2 -1 ' 

(A9) 



ng } (t,M 2 ,M 



/2a 



AT 
6M 2 Ji/M 



m 2 a(M 2 +aM /4 ) t( a M ,2 -l) 
(^Q/ g M 2 (aM' 2 -l) g M 2 + M> 2 . 



a 



(aM 1 



/2\2 



m 2 c a(M 2 + M ) 



, 2 m 2 c a 2 M'\M 2 + M' 2 ) tM 2 (aM rz - 1 



/2\ 



r/2 



- 2M Z - M - 



M 2 (aM' 



+ 



M 2 (aM' z - 

m 2 c a 2 M' 2 (M 2 + M' 2 )(M 2 + 2M' 2 - aM' 4 ) tM' 2 (aM 2 + 1 



(M 2 + M' ) 



+ 



M 2 (aM' 2 - l) 2 



(M 2 + M' ) 



+ 



+ 



2aM 2 M' 2 + M' 2 + aM' 4 
aM' 2 - 1 



(A10) 
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